Astronomy & Astrophysics manuscript no. 15346 


©ESO 2010 


November 5, 2010 





Modeling of the Vela complex including the Vela supernova 
remnant, the binary system y 2 Velorum, and the Gum nebula 

I. Sushch 1 ' 2 , B. Hnatyk 3 , and A. Neronov 4 



1 Humboldt Universitat zu Berlin, Institut fur Physik, Berlin, Germany 

2 National Taras Shevchenko University of Kyiv, Department of Physics, Kyiv, Ukraine 

3 National Taras Shevchenko University of Kyiv, Astronomical Observatory, Kyiv, Ukraine 

4 ISDC, Versoix, Switzerland 



Received 07 July 2010; accepted 18 October 2010 



O 

(N 

> 
O 

< 

O 

6 

-I— > 

c§ 1 1. Introduction 



ABSTRACT 



We study the geometry and dynamics of the Vela complex including the Vela supernova remnant (SNR), the binary system y 2 Velorum 
and the Gum nebula. We show that the Vela SNR belongs to a subclass of non-Sedov adiabatic remnants in a cloudy interstellar 
medium (ISM), the dynamics of which is determined by the heating and evaporation of ISM clouds. We explain observable charac- 
teristics of the Vela SNR with a SN explosion with energy 1.4 x 10 50 ergs near the step-like boundary of the ISM with low intercloud 
densities (~ 10~ 3 cm" 3 ) and with a volume- averaged density of clouds evaporated by shock in the north-east (NE) part about four 
times higher than the one in the south-west (SW) part. The observed asymmetry between the NE and SW parts of the Vela SNR 
could be explained by the presence of a stellar wind bubble (SWB) blown by the nearest-to-the Earth Wolf-Rayet (WR) star in the 
y 2 Velorum system. We show that the size and kinematics of y 2 Velorum SWB agree with predictions of numerical calculations for 
the evolution of the SWB of M ini = 35M Q star. The low initial mass of the WR star in y 2 Velorum implies that the luminosity of the 
nuclear line of 26 Al, produced by y 2 Velorum, is below the sensitivity of existing gamma-ray telescopes. 
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The Vela complex is one of the most interesting regions in the 
galactic plane. The observable flux from the Vela region ranges 
from radio to TeV energies. It consists of many objects, includ- 
ing the Gum nebula, the Vela supernova remnant (SNR) and Vela 
Jr. SNR (SNR RX J0852.0-4622 superposed on Vela), the binary 
system y 2 Velorum, the IRAS Vela Shell, an OB-association, £ 
Pup, etc. Some of them are shown in Fig. [2] Because of the dis- 
tances to these objects it is possible that some of them can in- 
tersect not only in projection, but also physically. We propose a 
scenario of joint evolutionary interaction of the Vela SNR, the 
binary system y 2 Velorum, the IRAS Vela Shell, the Vela OB2- 
association, and the Gum nebula. 



2. A hydrodynamical model of the Vela SNR 
evolution 

2.1. Vela SNR properties and peculiarities. NE-SW 
asymmetry. 

The Vela SNR is one of the closest supernova remnants to us. 
Different estimates of the distance to the Vela SNR su ffer from a 
large uncertainty: fro m 250 + 30 pc dCha et all 1 1999h to 350 pc 
dDubner et al.l(ll998h and references therein). Hubble Space tele- 
scope parallax observations of the Vela pulsar give the distance 
to the pulsar of D Ve ia = 294+™ pc dCaraveo et alll200lh. and the 
best e stimate is from the VLBI parallax measure dDodson et all 
l2003h : 



Similarly uncertain are the estimates of the Vela SNR a ge, 
which range from a few thousand years (|StothersL Il980l) to 
f SNR 2.9 x 10 4 vr (lAschenbach et all [1995^ 7 The most com- 
monly cited estimate is ?snr = W sar - 1-14 x 10 4 yr, w here 
fpuisar is the age of the Vela pulsar dReichlev et all Il970h . The 
total 0.1 - 2.4 keV X-ray luminosity from Vela SNR in erg/s is 
dLu & Aschenbachll2000l) 



L x , ot = 3.0xl0 35 



A 



290pc 



erg/s. 



(2) 



The main peculiarity of Vela SNR is the difference in the 
X-ray brightness and radius of its south-west (SW) and north- 
east (NE) parts. The ROSAT All-Sky Survey image (Fig. Q} 
of the Vela SNR reveals a shell with a diameter of about 8° 
dAschenbach et all 1 19951) . which implies a mean linear diameter 
of 



rfvela ~ 40 



D 



Vela 



290 pc 



pc. 



(3) 



The SW part of the shell appears to have a radius larger by a 
factor of 



/?sw - 1.3-Rj, 



(4) 



Dveia = 287^ pc. 



(1) 



than the NE part. 

Apart from the differe nce of the radii, the spatia lly-resolved 
spectroscopic analysis bv iLu & Aschenbac hi d2000h shows that 
the SW part of the shell appears to be hotter than the NE one. 
The shell is bright only on the NE side, while the SW side ap- 
pears to be dim and is apparently more extended (see Fig. [TJ. 
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The bounda ry between the bright an d the dim part of shell is 
quite sharp dLu & Aschenbachll2000l) . The change in the prop- 
erties of the shell at different sides indicates that the characteris- 
tics of the ISM in the NE part differ from those in the SW part. 
The sharpness of the boundary shows that the change in the ISM 
pro perties is abrupt rather tha n gradual. 

(Lu & Aschenbachl d2000t) estimate the contrast change of the 
two faint regions, each 1.5° x 1.5°, in the NE and SW parts, to 
be a factor of 1 1 in brightness and ~ 6 in emission measure. 
Assuming the same estimate for the emission measure contrast 
change for the entire NE and SW parts of the SNR and taking 
into account the difference of the radii of the two parts of the 
SNR, one can find that the luminosity of the NE part L x> ne is a 
factor of ~ 3.7 higher than the luminosity L Xt sw of the SW part. 

Apart from the difference in the overall luminosity and tem- 
perature, the characteristics of the X-ray emission from the NE 
and SW parts of the SNR also reveal a difference in the col- 
umn density of the absor bing material along the line of sight 
dLu & Aschenbachl l2000h . The absorption column densities Nu 



range from 5.0 x 10 cm in the NE part to 6.0 x 10 cm"" in 
the SW part. 

Finally, the Vela SNR is peculiar in still another aspect: the 
main shock of the SNR is not observed. Instead, the bulk of the 
X-ray emission is distributed all over the SNR volume. Such an 
observational appearance can be caused by the SNR expand- 
ing into a highly inhomogeneous ("cloudy") medium. In this 
case, the main shock advances through a low-density interstellar 
medium (ISM), leaving behind denser clouds, which are subse- 
quently heated and partially evaporated by thermal conductivity 
and transmitted shocks. This results in the appearance of a dis- 
tributed emission throughout the entire volume of the SNR, in- 
stead of from a thin shell at the interface of the main shock with 
the low-density ISM. The main X-ray emitters in the remnant 
are the two (cool and hot) phases (components) of heated cloud 
matter dLu & Aschenbac hL hoOOtlMiceli et al.Ll2005ll2006l) . 



2.2. A hydrodynamical model of the Vela SNR 

Despite extensive investigations, there is no common agreement 
about the evolutionary status of the Vela SNR. In early studies, a 
distance of 500 pc was assumed, which corresponds to a rather 
large size of 70 pc in diameter. Together with this assumption, 
the absence of the X-ray limb-brightening effect and of the clear 
signature of the high-speed shock in the SNR suggested a ra- 
diative stage of the Vela SNR evolution. The appearance of fil- 
ament structures in optic and radio waves as well as 100 km/s 
absorption lines in the spectra of background s t ars ad ditionally 
supported this scenario (see, e.g. iGvaramadzd (1 19991) and ref- 
erences therein). But recent observational data strongly support 
the hypothesis of the Vela SNR being in the adiabatic stage. New 
results for the distance place the Vela SNR at D « 290 pc, imply- 
ing that the mean radius of the Vela SNR is R « 20 pc, and the 
dynamics of bullets (ejecta fragments) outside the SNR bound- 
ary suggest the expansion of the SNR in the low-density ISM. 
The presence of the Vela pulsar indicates that the Vela SN pro- 
genitor was of Mm — 25 M Q . For M- m j = 1 1 - 25 M presuper- 
nova masses are Mf „ - 10.6- 16.6 Mr and m a sses of the ejecta 
M ei = 9 -H 15 Mr,. dLimongi & Chieffil d2006h . lEldridge & Toul 
d2004l) . iKasen & Wooslev et alJ d2009h V Therefore, the interac- 
tion (merger) of the massive ejecta with a velocity of over 1000 
km/s with the RSG shell occurs in the adiabatic regime, and the 
Vela SNR with a radius of 20 pc and an age of about 10000 
years should be in the adiabatic stage of evolution with char- 
acteristic velocities of about mRsNu/tsNR ~ 1000 km/s (m < 1 




09h15min 08 h 50 min 08 h 25 min 08 h 00 min 
RA (2000) 



Fig. 1 . ROSAT All-Sky Survey image (0. 1-2.4 KeV) of the Vela 
SNR dAschenbach et all Il995h . A-F are extended features out- 
side the boundary of the remnant ("bullets"). Light blue to white 
contras t represents a contrast in surface brightness of a factor 
of 500 dAschenbach et al.L[l995l) . Blue curves show the NE and 
SW hemispheres of the Vela SNR. The yellow curve shows the 
contour of the SWB of y 2 Velorum. 



in the free expansion case and m = 0,4 in the adiabatic one) 
without forming a thin dense radiative shell. Another observa- 
tional confirmation of the absence of the 100 km/s shell follows 
from the studies of absorption lines in the spectra of background 
stars dCha&Sembachil2000T) . where all stars in the Vela SNR 
direction at distances smaller than 350 pc do not show evidences 
of a 100 km/s absorption line. The maximum broadening corre- 
sponds to < 50 km/s. Only stars with distances exceeding 500 pc 
show 100 km/s features. Radio and optical shells (shell-like fila- 
ments) can be naturally explained by the emission of filamentary 
structures, exited by the SNR shock (akin to the Cyg Loop case, 
where filaments coexist with the main fast adiabatic shock). 

The hydrodynamical model of the Vela SNR evolution can - 
not be described directly by the Sedov solution dSedovL Il959l) . 
because the expansion proceeds in a cloudy rather than homo- 
geneous ISM. Instead, to describe the Vela SNR evolution one 
can use the White & Long solution, which describes the evo- 
lution of a supernova remnant expanding into a cloudy ISM 
dWhite&Longlll991l) . 

As explained above, the observed asymmetry of the Vela 
SNR is most probably due to the difference in the properties of 
the ISM on different sides of the remnant, in particular, by the 
different densities of clouds or different density contrast between 
the clouds and the IC medium. Below we assume that before the 
SN explosion the average number density (concentration) and, 
therefore, volume filling factor of clouds in NE part was larger 
than that in the SW part. This naturally explains the smaller ra- 
dius and the lower temperature of the NE part of the remnant. 
Once more, interaction of shock wave with clouds results in nu- 
merous filamentary structures of disrupted cloud material, visi- 
ble in radio band, and the radio image of the Vela SNR really 
shows more numerous shell structures in the NE part than in 
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SW one (iBock et all Il998l) . The NE-SW asymmetry proposed 
here can also explain a pulsar wind nebula (PWN) displacement 
to the SW according to the pulsar position. Namely, the reverse 
sho ck from more dense NE part reaches and de stroys PWN ear- 
lier dBlondin et all 1200 U lLaMassa et alll2008l) 

The conjecture that the densities on the two sides of 
the SNR are different is supported also by the dynam- 
ics of the shrapnel (bullets) - high-velocity clumps of 
the SN ejecta w it h an overabundan c e of heavy elements 



dAschenbach et all [1991 iMicelietall 120081: lLaMassa etafl 
120081; lYamaguchi & Katsudal 120091) . X-ray spectra of shrapnel 
argue a low density of the ISM around the Vela SNR, whereas a 
difference in the distances traveled by the protruding shrapnel on 
the NE (shrapnel pieces A-D/D') and on SW (shrapnel pieces E 
and F) sides of the SNR, as it is illustrated for the shrapnel pieces 
D and E shown in Fig.[T] is consistent with a corresponding dif- 
ference in mean ambient density, assuming that each shrapnel 
piece has approximately the same density and initial velocity. 

The solution for the shock radius r, of the remna nt as a func- 
tion of the age t has the form (I White & Lond.ll 99 ll) 



25(y + l)KE 



I6np ic 



1/5 



,2/5 



(5) 



where E is the explosion energy, p K is the density of the inter- 
cloud medium and y — 5/3 is the adiabatic index. The solution 
(0 differs from the standard Sedov solution without evaporating 
clouds by a different choice of a phenomenological constant K. 

The constant K depends on two parameters: on the ratio C 
of the fraction of the ISM mass initially contained in clouds and 
evaporated behind the shock wave to the ISM mass in the inter- 
cloud medium 



C 



(Pc) 
Pic ' 



where (p c ) is the volume-averaged density of the clouds, and on 
the ratio of the cloud evaporation time scale f ev to the SNR age t 



T = f ev /f. 



(7) 



To express the dependence of K on these two parameters we can 
use a simple analytical approximation 



K 
K~s 



(8) 



and K/K s ~ CT 1 for C « 1, t » 1, where K s = 1.528 is the 
value of constant K for the Sedov solution (when C = 0). From 
eq. ([5} follows the value of the shock velocity 



2 r s 
57 



(y + l)KE 



1 1/2 



which allows to calculate the temperature at the shock 



2(y- l)pm H 
(y+1) 2 k 



(9) 



(10) 



(here p. = 16/27 is the average mass per particle in hydro- 
gen mass units m# for helium/hydrogen abunda nce ratio »He - 
(1/12) »h and k is the Boltzmann constant) for the lWhite & Long! 
dl99lf) solution. 

The dynamics of expansion of the SNR is fully characterized 
by a set of four parameters: the explosion energy E, the preshock 
intercloud ISM density pj c , the cloud/intercloud density ratio C 
and the evaporation time/SNR age ratio t. The values of these 



parameters can be derived from the set of the observed charac- 
teristics of Vela SNR, such as the shock radii, characteristics of 
the X-ray radiation for the both NE and SW parts, etc. 

For the NE and SW shock radii /?ne — 18 pc and /?sw - 
23 pc and SNR age f Ve ia - 1-14 x 10 4 yr from eq. (0 we obtain 
corresponding shock velocities 



V NE = 0.4/? NE /f V ela ~ 6 -° X 107 Cm / S ' 

y sw = 0.4fl sw /fveia ~ 7.7 x 10 7 cm/s, 



(11) 

(12) 



niV£ 

s 



4.8 x 10 6 K, T s s w 



and from eq.dTOb shock temperatures T$ 

7.8 X 10 6 K. Average (emission measure weight ed) temperature 
of plas ma inside the SNR is about twice as high dWhite & Lond 
(fl99~flV ): 



< t£ >« 

< Tg >= 



l.9T" E ^9x 10°K, 



l.9T§ W : 



1.5 x 10 7 K. 



(13) 
(14) 



According to the results of & Asch enbac 2 (Eoooh X- 
ray radiating plasma in both parts of the Vela SNR consists of 
two phases: a hot one (T ~ 0.5 — 1.2 keV) and a cool one 
(T ~ 0.09 - 0.25 keV), and the cool one dominates the X-ray 
luminosity of the SNR. The temperatures found above corre- 
spond to the hot component in the two-temperature Raymon d- 
Smith thermal plasma model used in lLu & Aschenbachl (120001) '). 
The hot evaporated gas component with the volume filling factor 
fhot — 1 dominates the shock dynamics, while the cool one with 
fcooi = 1 - fhot ^ 1 dominates in X-ray radiation. The role of 
the initial intercloud interstellar gas is negligible in both shock 
dynamics and X-ray radiation. It means that 



(6) i 



hot 



» 1, 



where C hot =< p cJw , > /p u 
the approximation K/Ks 
radius is reduced to 



(15) 

and hereafter we take r <K 1 and use 
In this case, eq.® for shock 



^hof 



'"Vela = 



25{y+\)K s E 



\6nm H n hl 



1/5 



,2/5 



(16) 



where tihot - pc,hotl m H is the nucleon number density of the 
hot component. As we can see, Eq.(IT6b is the Sedov solution, in 
which a mean density of intercloud plasma inside the remnant 
Pi c is replaced by the mean density of the evaporated clouds (of 
the dominant hot component p c j, ol in our case). 

The nucleon number density of the hot component «/ !0 , can 
be estimated from its X-ray radiation. The X-ray luminosity of 
the SNR is an integral over SNR volume *V for X-ray emissivity 
ex = n e riHAx(T) of plasma with temperature T, electron (hydro- 
gen) number density n e (nu) and cooling function Ax(T) 



-L- 



exd'V. 



(17) 



Within the model of IWhite & Lon el (119911) . L x is estimated in 
erg/s as 



1.7xlO J4 eA_ 22 



1 cm 



1 pc 



1 + 



C 



1 +T 



(18) 



where n lc = pic/mn is the intercloud ISM nucleon number den- 
sity (nucleon number density n = (4/3)«h )), A_22 is the cooling 
function in units 10~ 22 ergcm 3 s~', Q is a number on the order 
of one and depends on C and t (for C » 1 and t <k 1 Q =* 1). 
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Unfortunately, we do not know the X-ray luminosity of 
the hot component of Vela SNR with satisfactory accuracy. 
Therefore, we use the more accurate data about the emission 
measure of the hot plasma for the calculation of the nucleon 
numb er density of the hot compo nent. According to the observa- 
tion of lLu & Aschenbachl(l2000l) . the mean value of the emission 
measure A of (p P i xe i x <p P i xe i = 8.75' x 8.75' pixel 



A = 10" 



f n e 

J pixel 



n H d<VI4nD%^ 



(19) 



of the hot plasma in NE region is Av~ - (5 -f 7) X 10 cm s . 
Meanwhile, from eq. ( fT9l and conditions n e nn = 0.66« 2 , 
'V pixel = (VpixelDvdaj 2 x 2R SNR it follows that 

r ff r 1/2 

nf 1 ' 2 = 1.5A 1 ' 2 ^ ™- 3 



20 pc 

and, taking into account that f. 



(20) 



■NE 
hot 



nil ^4.0x 10" 2 cm" 3 . 



Now, from eq.([T6b, we can find nucleon number density of 
hot component in SW part of Vela SNR 



SW _ NE 
n hot ~ n 



'hoi 



a 1 x 10" 2 1 



(22) 



Table 1. Characteristics of Vela SNR 



Parameter 


NE 


SW 


E, [erg] 


0.14 x 10 M 


t, [years] 


11400 


«> LpcJ 


18 


23 


Vc Mem /si 


600 


770 


T s , [K] 


4.8 x 10 s 


7.8 x 10 6 


nhot, [cm -3 ] 


0.04 


0.01 


fhot 


0.93 


0.91 


T llo „ [K] 


9x 10 6 


1.5 x 10 7 


n cooh [cm" 3 ] 


0.38 


0.10 


fcool 


0.07 


0.09 


T C ooi, [K] 


1 x 10 6 


1.7 x 10" 



where we use the approximation validated by IWhite & Long! 
( fl99lh A_ 22 (r) = 1. Similarly, 



L lot - 0.4 xl0 35 A_ 22 ergs" 1 , 



(21) L, 



cool 



0.6 x 10 35 A_22 ergs" 1 , 



L S Z - 0.06 x 10 35 A_ 22 ergs" 1 . 



(27) 
(28) 
(29) 



The total model luminosity Lx,tot - 3.2x 10 35 A_22 ergs 1 is close 
to the observed one of eq. (O 



and the energy of the Vela SNR explosion 
l6mnH 



E = 



25( r +l)^4 



1.4 x10 s0 



erg. 



(23) 



The cool component of the X-ray radiating plasma in the 
Vela SNR does not influence the shock dynamics, but dominates 
in the X-ray luminosity. The mean value of the emission measure 
A and the temperature of the cool plasma in the NE region is 
A N J ol - (3 t 5) x 10~ 3 cm" 5 , T NE , - 0.1 x 10 7 K and in the 



SW region A™, 



cool 
1-4, 



(4 -r 6 ) x 10 4 cm- 5 , T sw , - 0.17 x 10 7 K 

■ 1- VVIM 1 | ^ ■ y ' COOl 

(iLu & Aschenbachi 120001) . From eq.© for the NE part of the 
Vela SNR follows 



("coolj c00 i) 



1-6(A^,) 



l ' 2 cm" 3 



1.0 x 10-' cm" 3 . 



(24) 



The filling factors of hot and cool plasma in the NE part can 
be estimated from the equality of the pressure P oc nT of both 
components: (n coo iT coo i) NE - (n hot T hot ) NE 



sNE 
J cool 

fNE 
J hot 



fNE 
J cool 

^ ^ cool 



T coo i( n coolf r l„i) 



1 coor 



1 hot WhotJhot > 



- 8 x 10- 



(25) 



or f* o E ol - 0.07, ff E t - 0.93. And, finally, from eq.® follow 



nucleon number densities n NE , 

cool 



0.38 cm 



3 <t 



0.04 cm- 3 . 



Similarly, for the SW part of the Vela SNR (n coo if coo 
3.1 x 10- 2 cm- 3 , fJ ol - 0.09, fj - 0.91, and n s J ol - 
0.10 cm" 3 , nH ^0.01 cm- 3 . 

Table 1 gives a summary of the parameters of the NE and 
SW parts of the SNR derived from the X-ray data. 

As a test for the self-consistency of our model we 
can calculate the predicted X-ray luminosities of different 
parts/components of Vela SNR in erg/s, using eq.dT7b: 



j NE _ ^ r>3 
L cool - 3 nR NE 



f E ol (n e n H A x (T))Z 



- 2.2 x 10 35 A_ 



22. 



(26) 



3. Interaction of the Vela SNR and the y 2 Velorum 
stellar wind bubble 

Estimates of the physical parameters of the Vela SNR derived in 
the previous section show that it is evolving in the inhomoge- 
neous medium whose density changes in a step-like manner. In 
this section we explore the possibility that this step-like change 
of the ISM parameters can be related to the presence of the 
boundary of the stellar wind bubble (SWB) around a Wolf-Rayet 
(WR) star in the y 2 Velorum system, which is situated in the 
vicinity of the Vela SNR. 

3. 1 . Physical characteristics of the y 1 Velorum binary 
system. 

y 1 Velorum is a WC8+08-8.5III binary system. The WR com- 
ponent in this binary system is the closest to Earth WR star 
(WR11). The estimates of the distance to y 2 Velorum have 
evolved over the years and are still controversial. An early 
estimate o f the distance based on HIPPARCOS parallax was 
258^ j pc dSchaerer et all 1997[). It was rece ntly revised by sev- 
eral scientific groups. iMillour et al.l (120071) give an in t erfero- 
metric estimate of the distance 368^ pc. iNorth et all (120071) 
estimate the distance to y 2 Velorum to be 336^ pc based on 
the orbital solution for the y 2 Velorum binary obtained from 
the interferometric data. Finally, a revision of t he analysis of 
HIPPARCOS data gives a distance of 334^° pc Ivan Leeuwenl 
d2007h . Therefore, in the following we take the distance to 
y 2 Velorum to be D y i Vet - 330 pc. 

The current mass estim ate of the WR star is Mwr = 9.0 + 
0.6M o dNorfh et al. 



.- 120071) The mass of the O star is M = 
(28.5 ± l.l)M n (INorth et al Ll2007b . The y 2 Velorum system is an 
important source from the viewpoint of nuclear gamma-ray as- 
tronomy. Because it is the nearest WR star, this is the only source 
that can potentially be detected as a point source of 1.8 MeV 
gamma-ray line emission from the radioactive 26 Al with current 
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generation instruments. Previous observations by COMPTEL 
put an upper limit at the level of 1.1 x 10 5 J cm _2 s _1 on the 
line flux from the source dOberlack et which is, appar- 

ently, below the typical predictions of the models for stars with 
the initial mass M m \ ~ 6OM . One should note, however, that 
modeling of 26 Al production suffers from uncertainties of the nu- 
clear reaction cross-se ctions, stellar parameters (such as rotation 
and m etallicity), etc. dLimongi & Chieffil l2006t iPalacios et al.L 
l2005h . 

The estimate of the initial mass of the WR star, Mj n ; ~ 
(57+15)M , is obtained from the evolutionary models of isolated 
rather than binary stars dSchaerer et al.L IT997h . It is possible that 
within a binary system, the mass transfer between the compan- 
ions can change the stellar structure of both components, so that, 
for example, stars with initial m asses as low as 2OM can be- 
come WR stars in binary systems (I Vanbeverenl 1 1 99 ll) . Modeling 
the binary evolution of y 2 Velorum leads to a lower limit on the 
initial mass of the WR star: MwR,i > 38M in ( Vanbeveren et al.L 
11998b and M WR ,i - 35M Q , M ,i '- 31.5M in kldridgdl2009l 



3.2. The stellar wind bubble around y 2 Velorum 

WR stars are expected to be surrounded by the multi-parsec scale 
bubbles blown by the strong stellar wind inside a photoionized 
HII region. The size of the HII region and bubble depends on 
the (time-dependent) joint action of the photon luminosity in the 
Lyman continuum L^c and the mass-loss rate of the star, M w , on 
its age, f, on the wind velocity v w and on the density of the ambi- 
ent medium into which the bubble expands, po- Qualitatively, the 
relation between the bubble radius 7?bub and the a bove parame- 
ters ca n be found using the analytical calculation of lWeaver et al.l 
d!977h 



Kb 



1 



ub 



1/5 



Po 1/5 ^ 



(30) 



where k is a numerical coefficient which, in the simple case 
considered by IWeaver et all d 1977b was k = [125/(154tt)] 1/5 , 
while in a realistic case it can be found from numerical model- 
ing based on a stellar evolution mode l. For example, a numerical 
model considered by lArfhurl d2007l) . which studied the evolu- 
tion of a star with initial mass 40M o in a medium of density 
no = 15 cm 3 , predicts the final radius of the shock in the am- 
bient interstellar medium R s ^ =* 36 pc at the end of WR phase 
before SN explosion. Using the Eq. ( 130b one can re-scale the nu- 
merical simulations for the particular values of ISM density to 
find that for the typical value no =0.1 cirT 3 in the ISM the ra- 
dius of the SWB can reach 100 pc. Meanwhile, as we will see 
later, y 1 Velorum system was born in a molecular cloud and its 
bubble should be considerably smaller. 

3.3. Overlap between the Vela SNR andy 2 Velorum bubble 

Comparing the distances of Vela SNR (~ 290pc) and y 2 Velorum 
system (~ 330pc) and taking into account the proximity of the 
two objects on the sky, one can notice that if the SWB around the 
y 2 Velorum system has indeed a radius of 30-70 pc (> 5° on the 
sky), the Vela SNR, which itself has a size of ~ 40 pc, is expected 
to physically intersect with the y 2 Velorum bubble. In view of 
this geometrical argument it is natural to ascribe the observed 
step-like change in the parameters of the ISM at the location of 
the Vela SNR to the boundary of the SWB of y 2 Velorum. 

The hypothesis of intersection between the Vela SNR and 
y 2 Velorum SWB is further supported by the simple geometri- 



cal form of the boundary between the bright and the dim part 
of the SNR shell. Indeed, the boundary roughly follows the con- 
tour of an ellipse whose major axis is perpendicular to the direc- 
tion from the center of Vela SNR toward the y 2 Velorum, so that 
the minor axis is aligned with the direction toward y 2 Velorum, 
see Fig. Q] The projected distance between the Vela SNR and 
y 2 Velorum is D ' = 5.2°. Assuming the distances Dv e i a — 290 pc 
and D y i Vel =* 330 pc one finds that the physical distance between 
the two objects is 

Rf-Vel = £>Vela-y?Vel - 44 pc, (31) 

which we adopt as an estimate for the radius of the SWB around 
y 2 Velorum. 

3.4. Estimate of the total mass of the stellar wind bubble of 
y 2 Velorum 

Within the geometrical model discussed above, the observed dif- 
ference in the absorption column density Nh between the NE 
and SW parts of the Vela SNR can be used to estimate the to- 
tal mass of the ISM swept up by the stellar wind of y 2 Velorum 
over the entire lifetime of the SWB. Taking the difference be- 
tween the measured Nh values in the NE part and SW parts 
dLu & Aschenbachll2000b 



AAfe = AW - Nh.ne - 5.5 x 10 20 cnT 2 

one can find the total mass of the y 2 Velorum SWB: 



M - AnR 2 ylvd kN H X -m H « 1.3 X 10 s M Q 



R 



y 2 Vel 



44 pc 



(32) 



(33) 



Assuming that the bubble has expanded into a homogeneous 
ISM over the entire expansion history, one would estimate the 
ISM density around y 2 Velorum as 



M 



«ISM 



(4/3)n& 2Ve[ m H 



12 



R 



y 2 Vel 



44 pc 



cm 



(34) 



One could notice that this estimate of the density of the ISM 
is much higher than the estimates of both the intercloud and of 
the volume-averaged cloud density around Vela SNR (Table 1). 
We come back to the discussion of the origin of this discrepancy 
below. 



4. The y 2 Velorum system and the IRAS Vela shell 

The angular size of the y 2 Velorum SWB, found from our inter- 
pretation of the asymmetry of the Vela SNR, i.e., from the ratio 
of Ryiy e i — 44 pc and the distance 330 pc as st 7.6° coincides 
with the angular size of a large circular arc like structure, cen- 
tered on (l,b) = (263°, -7°) with the radius R lvs - 7.5°, which 
is visi ble in the inf rared band, known as the "IRAS Vela shell" 
(IVS) dSahuLll992l) . This structure surrounds the Vela OB2 asso- 
ciation including y 2 Velorum and f Puppis. A recent study of the 
spatial distribution of the neutral hydr ogen and radio conti nuum 
emission at 1420 MHz of the IVS bv ITestori etaD d2006l) pro- 
vides a new estimate of the coordinates of the centroid of IVS 
from the observed IR emission (l,b) = (259.9°, -8.3°) and the 
radius of the neutral hydrogen shell in the SW se ctor between 
positio n angles ~ 162° and ~ 265° Ri\s - 5.7°. ITestori et al.l 
d2006h also estimate the mass of ionized and atomic components 
of the shell, assuming a distance of 400 pc: M m - 9.1 x 1O 4 M 
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(or Mivs - 6.0 x 1O 4 M for our distance 330 pc). It is ex- 
pected th at the amount of molecular g as in the IVS is about 
~ 10 5 MM( lRaiagopal & Srinivasanl ( ll998l) ). 

A possible interpretation of the IVS as a boundary of the 
SWB of y 2 V elorum was discussed bv lOberlack et alJd2000t) and 
iTestori et alj d2006l) . 

Adopting this interpretation, one can find that the estimates 
of the parameters of the SWB of y 2 Velorum derived above from 
the analysis of X-ray data on the Vela SNR, agree well with those 
found from the analysis of the IVS data. In particular, adopting 
an estimate of the distance toward y 2 Velorum for the center of 
the IVS, one finds the radius 



R lvs ~ 44 



^Velorum 



pc 



330 pc 

and the total mass of the shell 
M IVS ~ 1 x 10 s M 



(35) 



(36) 



comparable to the estimates found in Eqs. OTI ) and d33b , respec- 
tively. 

4.1. Implications for the models of the Wolf-Rayet star 

The initial mass of the WR11 star M; n j is equal the sum of the 
current mass of the star Mwr = 9.0 M and the mass which was 
blown by the wind: 



M W r + AMmss + AM (RSG s/lbvs) + (m W r) ?wr, (37) 



where Mwr is the mass loss rate and twR is the time inter- 
val of the WR wind. The sum in the above equation includes 
mass loss via different types of winds ejected by the star at 
different stages of stellar evolution: the main sequence stage 
(MSS) wind, the red supergiant stage (RSGS) (for stars with 
initial mass M m \ < 4OM ) or the luminous blue variable stage 
(LBVS) (for Mi n i > 40M o ) w ind, and the continu i ng WR wind. 
Numerical calculations from Frever et al. ( 20031). Frever et al 



2006J) . Ivan Marie et all (12005 ). lArthud (l2007h . Ivan Marie et al 
(120071) . iPerez-Rendon et all 00091) show that the MSS and 



RSGS/LBVS dominate the mass loss with typical values of mass 
loss before WR stage (21 4- 26) M for M ini ~ (30 -H 60) M Q . 
Meanwhile, the WR stage dominates in the kinetic energy of the 
wind, injected into the wind bubble. 

For the y 2 Velorum binary one can calculate the total stellar 
wind mass loss AM 29M as the difference between the ini- 
tial and contemporary masses of stars in binary systems, which 
are M W R jini =; 35M Q , M ,ini - 31.5M and M W R _^9.OM , 
M - 28.5M respectively, according to lEldridgd (120091) . This 
includes the mass-loss rate from the slow red supergiant stage 
wind AMrsgs - 19M dFrever et al.l (120061) 1. Therefore, the re- 
maining 1OM of the hot intercloud gas in the y 1 Velorum SWB 
correspond to fast MSS and WR winds. Clumps, created by the 
interaction of the fast WR wind and the slow RSG wind, together 
with ISM clouds (Gum nebula interior, see below), survived the 
passage through the expanded IVS, are the main sources of X- 
ray emitting plasma in SW part of Vela SNR. As follows from 
Table 1, the lower limit on the total mass of gas in clumps and 
clouds inside the SWB is 



SWB 



M\ 



hot 



hoi fSW , cool rS W 
n SWB^hot + n SWB^cool 



256M . 



(38) 



Let us estimate the energy budget of IVS. Assuming that the 
IVS is a boundary of the SWB of y 1 Velorum we can estimate its 



kinetic energy, assum ing an expansion velocity V exp - 1 3 kms 
(ITestori et alj (120061) 1 



7 ivs 



^~M IVS V 2 -2xl0 50 erg. 



(39) 



Thermal energy inside of the SWB of y 2 Velorum is limited 
by conditions of strong Vela SNR shock and the dominance of 
the mean density of evaporated clouds in postshock plasma. 

Strong shock condition means that the ratio of pressures in- 
side the remnant and bubble and equivalently the ratio of energy 
densities should be high ER S 3 W /E^ VB R-^ S » 1 or 



E™ B « E 



R 



IVS 



R, 



1 x 10 51 ere 



(40) 



The dominance of evaporated cloud material in the postshock 
region Cf w = n sy f /n SWB » 1 means that n SWB < O.lnf w f and 

& hot c.,hot' ic ic hot 

we take hereafter n SWB = 10" 3 cnr 3 (i.e. C = 10) and T SWB = 

IC v 7 IC 

3 x 10 6 K as reasonable parameters of the hot (intercloud) gas 
inside the SWB enclosed in the IVS. The total mass and thermal 
energy of this hot gas inside y 2 Velorum is M S ^ B 1OM and 
E s h ™ B ^ 1 x 10 49 ergs. Meanwhile, the thermal energy of clumps 
and clouds with a reasonable value of temperature T c i < 10 4 K 
is only E s f B < 1 x 10 48 ergs. 

To summarize, in our hydrodynamical model of IVS the total 
(thermal and kinetic) energy of the y 2 Velorum SWB is 



p _ plVS , pSWB pSWB „ r. in5() 

Eswb - E kin + E hot +E cl ^2X10 erg 



(41) 



with evident dominance of kinetic energy of the massive (~ 
1O 5 M ) IVS in the total sum. Inside the SWB we predict about 
1OM of hot intercloud gas and about 26OM of immersed 
warm/cool clumps/clouds (in case of their complete evaporation 
inside Vela SNR). It is worth noting that temperature and density 
of plasma inside IVS are not well restricted by observations, and 
here we use approximate values, which are consistent with the 
adopted limits. 

We can compare these estimates with the results of nu- 
merical simulations of stellar wind bubbles around WR stars. 
Numerical simulations of the evolut ion of a s t ar wit h an initial 
mass M ini = 3 5M , as proposed b v lETdridgd (l2009l) for WR1 1. 
were made by iFrever et"al] d2006l) for an environment with a 
density of no = 20 crrT 3 and a temperature of To = 200 K. They 
show that at the end of the calculations (before the SN explo- 
sion) the hot gas bubble has a mean radius of 34 pc and shell- 
like HII and HI regions of the swept up ambient gas extend out 
to a distance of 43 - 44 pc, the total mass is 1 .5 x 1O 5 M , the ki- 
netic energy is 4.9 x 10 49 erg, the thermal energy of the hot gas is 
1 . 1 x 10 50 erg, of the warm gas 4.3 x 10 49 erg, i.e., the radius, total 
mass and total energy (kinetic and thermal, 2 x 10 50 erg) are sur- 
prisingly close to our estimate for y 2 Velorum. Nevertheless, two 
important differences should be clarified for y 2 Velorum namely, 
the kinetic energy dominance and the low density of the ISM. 



4.2. Density of the interstellar medium and interaction with 
the Gum nebula 

The three estimates for the "typical" ISM density derived above 
provide widely different values. The estimate based on the dy- 
namics of expansion of Vela SNR suggests a low value for the 
ISM density «ism ^ 0.01 cm 4 , while the estimate based on the 
total mass of the SWB around y 2 Velorum («ism ~ 16 cirT 3 ) and 
the estimate based on the dynamics of expansion of the SWB 
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Gum nebula as a fossil Stromgren sphere of the Vela X super- 
nova and estimated the distance to nebula to be 



270 



260 250 240 

Galactic Longitude (deg.) 



Fig. 2. Locations of the Vela SNR (Vela pulsar is shown as a 
cross), y 2 Velorum (shown as a circle), IRAS Vela Shell (IVS) 
bubble and Gum nebula (center is shown as a square) in Galactic 
coordinate system. 



("ism ~ 20 cirT 3 ) suggest a much higher density. This points to 
the fact that the distribution of the ISM in the direction of the 
Vela region, in the distance range ~ 300 - 400 pc, is highly in- 
homogeneous. This is, in principle, not surprising, because the 
region is known to contain several stellar formations with differ- 
ent properties. 

First, y 2 Velorum belongs to the yVel asso ciation, which 
is a s ubcluster of the OB-association Vela OB2 dJeffries et al.L 
120081) . The density of the ISM in the OB association soon 
after y 2 Velorum' s formation is expected to have been much 
higher because of the presence of a parent molecular cloud. 
The initial expansion of the y 2 Velorum bubble into a dense 
(lO 2 - 10 3 cm 4 ) molecular cloud (the progenitor of Vela OB2) 
can explain the large mass of the swept up ISM in a shell around 
the SWB. In this scenario, the stellar wind and the radiation of 
y 2 Velorum destroyed the parent molecular cloud and swept up 
its gas. For some time the stellar wind of y 2 Velorum was practi- 
cally trapped inside a ~ 1O 5 M cloud, and only when the energy 
accumulated in the SWB and HII region exceeded the gravita- 
tionally bound energy of the cloud (~ 10 50 ergs for R c i = 10 pc 
cloud), the dense shell of swept-up cloud material began to be 
accelerated by the thermal pressure of the SWB and the HII 
region gas without considerable additional mass loading and 
counter pressure of hot rarefied gas of the ISM. At this evolu- 
tionary stage, the thermal energy of the system converts into the 
kinetic energy of the shell, resulting in the atypical dominance 
of the kinetic energy of the IVS in the total energy balance of the 
y 2 Velorum SWB/HII region. 

Next, both the Vela SNR and the y 2 Velorum SWB could be 
interacting with a still larger scale SNR, known as the Gum neb- 
ula. This nebula is a very large region of the ionized gas about 
36° in diame ter, centered app roximately at (/, b) = (258°, -2°) 
(lGumlll952h . shown in Fig.l2l lBrandt et all (1197 lh explained the 



£>Gum =400±60pc. 



(42) 



iRevnoldsl (1 19761) suggested that the Gum nebula is an 1 
Myr old SNR, which is now heated and ionized by the two 
very hot stars ze t a Pup pis and gamma Velorum within it. Later 
IWoermann et alJ (1200 lh showed that the Gum nebula can possi- 
bly be a SNR of the zeta Puppis companion. They showed that 
runaway O-star zeta Puppis was within < 0.5 deg of the expan- 
sion center of the Gum nebula about 1 .5 Myr ago, which is ev- 
idence of the relation between the SN explosion of the binary 
companion and the Gum nebula expansion. Assuming a distance 
of 400 pc the radius of Gum nebula is about 



£>Gu 



400 pc 



pc. 



(43) 



This means that both the Vela SNR and y 2 Velorum are situated 
inside the cavity formed by the expansion of the SNR associated 
to the Gum nebula. 

The Gum nebula as a very old (0.9 - 2.0 Myr) SNR should 
be at the late radiati ve stage of evolutio n, which can be modeled 
with the solution of ICioffi et al.l (1 19881) for the shock dynamics. 
Assuming the the current radius 124 pc, age ~ 1.5 Myr and the 
explosion energy of 10 51 £5i erg, we can find out the density of 
the ISM in which the Gum nebula expands: 



"ism = 0.07? 



t 


6 


^Gum 


1.5 Myr 




400 pc 



£ 5 7 jcnT 3 , 



(44) 



where £ is the metallicity factor, equal to 1 for solar abundances. 

The estimate of the average density of the ISM around Gum 
nebula shows that actually the nebula expands into the medium 
with a much lower density than the one suggested by the estimate 
of the total mass of the SWB around y 2 Velorum. This supports 
the hypothesis that the enhancement of the density of the ISM 
was locally present around the y 2 Velorum system at moment of 
its birth, most probably because of an OB association. 

5. Conclusions 

We developed a model for the interaction of the Vela SNR 
and the y 2 Velorum SWB, which explains the observed NE/SW 
asymmetry of the Vela SNR. 

Adopting a model of the expansion of the Vela SNR into a 
"cloudy" ISM, we showed that the volume-averaged density of 
the shock-evaporated clouds in the NE part of the SNR has to be 
about four times higher than in the SW part. We noticed that a 
plausible explanation for the observed density contrast is that the 
Vela SNR exploded at the boundary of the SWB around a nearby 
Wolf-Rayet star in the y 2 Velorum system, which is situated at 
approximately the same distance as the Vela SNR. 

Within our model of the interaction of the Vela SNR and the 
y 2 Velorum SWB, the difference of the spectral characteristics 
of the X-ray emission from the NE and SW parts of the remnant 
can be used for an estimate of the parameters of the y 2 Velorum 
bubble. We showed that the measurement of the change of the 
column density of the neutral hydrogen gives an estimate of the 
total mass of the SWB, ~ 1 x 1O 5 M . 

On the basis of modeling the dynamics of expansion of the 
bubble around y 2 Velorum, we confirmed the initial m ass of the 
Wolf-R ayet star in the y 2 Velorum system suggested bv lEldridgel 
(120091) to be = 35M . This estimate is lower than the previous 
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estimates used for the derivations of predictions of the flux of a 
y-ray spectral line at 1.8 MeV, expected from the decays of 26 Al 
in this source. Taking into account the revised estimate of the 
initial mass of the Wolf-Rayet star, the 26 Al line flux flux from 
y 1 Velorum is expected to be much below the COMPTEL limit. 
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